arXiv:l502.04274v2 [quant-ph] 24Jul2015 


A Fundamental Form of the Schrodinger Equation 


Muhammad Adeel Ajaib^ 

Department of Physics and Astronomy, Ursinus College, Collegeville, PA 19f26 


Abstract 

We propose a first order equation from which the Schrodinger equation can be derived. 
Matrices that obey certain properties are introduced for this purpose. We start by con¬ 
structing the solutions of this equation in ID and solve the problem of electron scattering 
from a step potential. We show that the sum of the spin up and down, reflection and 
transmission coefficients, is equal to the quantum mechanical results for this problem. 
Furthermore, we present a 3D version of the equation which can be used to derive the 
Schrodinger equation in 3D. 
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1 Introduction 


The problem of particle scattering from, for instance, a potential barrier or step are important 
fundamental examples in quantum mechanics. The Schrodinger equation can be employed to 
solve these problems in the non-relativistic limit. The Schrodinger equation, however, does 
not consider the spin of the particle in such scattering problems. For relativistic particles, the 
Dirac equation [1] when applied to such problems leads to the Klein paradox [2]. This paradox 
arises when the Dirac equation is used to study problems such as the scattering of an electron 
from a potential step, particularly when the barrier height Vq > 2mc^. For this problem the 
reflected current turns out to be larger than the incident current. This paradox can be resolved 
by assuming the production of particle-antiparticle pairs at the barrier. Other methods have 
also been proposed in literature to resolve this paradox (see, for example, the list of references 

[ 3 ]) 

In this paper we introduce a non-relativistic hrst order equation which can be employed to 
include the spin of particles in scattering problems. We hrst show that the Schrodinger equation 
can be derived from this fundamental hrst order equation in ID. We derive the probability and 
current densities for this equation. Next, we construct the solutions to these equations and 
calculate the rehection and transmission coefficients for this process. We show that combining 
the results for spin up and down electrons leads to the quantum mechanical result for the 
rehection and transmission coefficients. Our analyses makes specihc predictions about the 
manner in which spin up electrons can scatter oh a step potential. 

The paper is organized as follows: In section 2 we present a fundamental form of the 
Schrodinger equation in ID and derive the probability and current densities. We further discuss 
the plane wave solutions to this equation. With the solutions constructed, we solve the potential 
step problem for E > Vq and E <Vo in ID in section 3. We introduce the 3D version of this 
equation and derive the probability current in section 4. Our conclusions are presented in 
section 5. 


2 An Equation Underlying The Schrodinger 
Equation in ID 

In this section we present an equation that, even in the non-relativistic limit, can be seen as 
an underlying equation to the Schrodinger equation. The Schrodinger equation can be derived 
from this equation similar to the manner in which the Klein Gordon equation can be obtained 
from the Dirac equation. We restrict ourselves to a single spatial dimension (say, z) in this and 
the following sections. In section 4 we will present a 3D version of the equation. 

We propose the following first order equation 

— idzip = {i'r]dt + ( 1 ) 
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Here rj and 77^ are matrices that obey certain relations described below. Using iteration the 
above eqnation leads to the following 

{—idzfi) = {irjdtft!) + im {rj, //'*'} dt'ip + (2) 

In order to obtain the Schrodinger eqnation the matrices rj and 77^ should satisfy the following 
properties 

772 = 0 (3) 

{rff = 0 (4) 

{>?,’;*} = 2 / ( 5 ) 


In addition, we require the 77 matrices to be symmetric. We hnd that the minimum dimension 
for the matrices 77 that are symmetric and satisfy the relations above is four. Moreover, these 
matrices were found as pairs of two, one being the conjugate of the other. Therefore, the 
matrices 77 are symmetric 4 x 4 nilpotent matrices and the only eigenvalue of these matrices is 0 . 
Furthermore, the trace and determinant of these matrices are also zero. The wave function 7/’ 
in equation ( 1 ) is therefore a column matrix with four components. There are several possible 
representations of the matrix 77. For instance, two possible representations of this matrix are 
as follows 
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where cXi are the Pauli matrices. The matrix 77 is symmetric (77^ = 77) but not Hermitian, 
therefore, 77I = 77* 7^ 77. For our analysis herein, we choose the representation in equation (6). 
The choice of this representation is primarily because the eigenstates corresponding to this 
representation are more relevant for the analysis in the following sections. 


Probability Current 

Multiplying equation ( 1 ) with 7/'I(77 + 77I) from the left we obtain the following equation 

— i^jJ\r| + ri^)dz'4> = {iijj^ri^ridt)'4> + ( 8 ) 

Similarly, taking the complex conjugate of equation ( 1 ) and multiplying by (77 + on the 
right we obtain the following equation 

+ 77^)'0 = rf)il) + m'lp^ 7 ]rj^ ip ( 9 ) 
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Subtracting (8) from ( 9 ) we obtain the continuity equation 

( 10 ) 

with the probability and current densities given as 

J = ^jJ^ {r] + r]^)%jj ( 11 ) 

p = (12) 

The probability and current densities are therefore Hermitian. For the representation of rj in 
(6), T] + = —i\/2'j2 and = I + i'js, where y, are the gamma matrices. 

2.1 Plane Wave Solutions 

We hrst consider motion along a single {z) direction and seek plane wave solutions of equation 
(1) of the form 


ij = u{p)e-^P-^ = 


( 13 ) 


where, E = pl/ 2 m. Therefore, in momentum space equation ( 1 ) becomes 

Pz = (Ep + mp"') '0 


( 14 ) 


The above equation is an eigenvalue equation with the momentum operator given as, P = 
Ep + p^m, which is not Hermitian. However, the eigenvalues of this operator are real. The 
eigenvectors of the matrix Ep + mp'^ consists of two states with eigenvalues of momentum 
+y/ 2 Em = Pz and two other with states —\/ 2 Em = —pz- This is consistent with the energy 
momentum relationship of particles in quantum mechanics, E = pi/ 2 m. Following are the 
eigenstates of the matrix Ep + ylm. 
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( 15 ) 


( 16 ) 


where a = 1 /{E + m) and Pz = \/ 2 Em. The eigenstates corresponds to a spin up particle 
with positive momentum and corresponds to spin down particle with the same momen¬ 
tum. The eigenstates and correspond to particles with negative momentum, i.e.. 
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Figure 1: Particle incident on a step potential. 

moving in the negative z direction. The eigenstates are related as = u^^\—Pz) and 

u^'^KPz) = u^'^K~Pz)- We will employ these eigenstates in the following section when we study 
the scattering problem. The normalization condition for the eigenstates are 

= 2 (17) 

^(l)t^(2) ^ Q ('Ig) 

and similarly for and Since the energy momentum relationship of these states does 
not follow the relativistic expression, it is not relevant to discuss how these states transform 
under Lorentz transformations. 


3 Potential Step problem 

The potential step problem is one of the foundational examples in quantum mechanics (Figure 
1). As described earlier, solving this problem with the Dirac equation leads to the Klein 
paradox. In this section we analyze this problem for electron scattering from a potential step in 
1 dimension using the solutions presented in previous sections. We will show that the resulting 
transmission and reflection coefficients for spin up and down particles are related to the ones 
in quantum mechanics. For regions I and II, equation (14) is given by 

Pi Ip ={Ep + mp^) Ip (/) (19) 

P2 i> ={{E - Vo)p + mp^) xp {II) (20) 
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3.1 Case I: £” > Vq 

We first consider the case of a spin up electron incident on a potential step with energy E > Vq. 
The incident and reflected electron waves in region I are given by 
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Similarly for region II we have 
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( 21 ) 


( 22 ) 


(23) 


where a = 1/{E + m), (3 = 1/{E — Vo + m)^ pi = V2Em and p 2 = ^/‘2{E — Vo)m. At = 0 
the continuity of the wave function implies 


%Iji{z = 0 ) + tjj'jiz = 0 ) = ilJn{z = 0 ) 

From this condition the coefficients are found to be 
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Figure 2: The plot shows the transmission and reflection coefficients for spin up and down 
electrons given in equations (A-1), (A-2), (A-3) and (A-4). The red lines show the coefficients 
of spin up electron and the blue lines correspond to spin down electron. The value of the 
potential step is chosen to be Vq =100 eV. For i? ~ Vq there is no transmission and the electron 
is completely reflected with its spin up. For E slightly greater than Vq, the transmission of spin 
up and down electron occurs with nearly equal probabilities. For E/Vq > 1.05 the transmitted 
electron is dominantly spin up. The lower panel shows the reflection coefficients. We can 
see that there is a very small probability that the reflected electron has its spin flipped. The 
reflection coefficient for spin down is small but non-zero. As the energy increases the barrier 
appears more and more transparent to the electron and the transmission coefficient of the spin 
up electron becomes large. For E ~ 3Vo the barrier is nearly transparent to the electron. 
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Figure 3: Plot shows the reflection coefficients (A-3) and (A-4) for Vq = 100 keV and Vq =1 
MeV. The red lines correspond to spin up electrons whereas the blue lines are for spin down 
electron. The plot of the transmission coefficients for these cases are similar to the one shown 
in Figure 2. We can observe that the reflection coefficient changes notably as the height of the 
barrier is increased. For Vq =100 keV, the probability of the spin up electron being reflected as 
a spin up or down electron is nearly equal. For Vq =1 MeV, the electron is dominantly reflected 
with its spin flipped. 
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The transmission and reflection coefficients Ti, T 2 , Ri and R 2 are given by 
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(30) 


The complete expression for these coefficients are presented in the appendix in equations (A-1), 
(A-2), (A-3) and (A-4). The sum of these coefficient is always equal to 1, i.e., 


(Ti + T2) + {Ri + R2) — 1 


(31) 


The quantum mechanical transmission and reflection coefficients are related to these coefficients 
as 
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(33) 


We can see that although the reflection and transmission coefficients (Ti, T 2 , -Ri and R 2 ) 
depend on the mass of the incident particle, the sum, which is the quantum mechanical result, 
is independent of it. The current densities for the incident, reflected and transmitted spin up 
and down particle are given by 


jj = = 

Jr = 

Jr = ^hiv + V^)^ri -■ 
Jt = + 

The conservation of probability density implies 


4|Appi 

E + m 


4pi|Bp 


E + m 
Api\BY 

2 

E + rn 
Ap2\C 

? 

E + m - 
4p2|C" 

-Vo 

2 

E + m - 

-K, 


(34) 

(35) 

(36) 

(37) 

(38) 


W! 4. hil ^ 

I Tt| I Tt| I Tt| I Tt| 

I'^/l r/l I'^/l l-Jil 


(39) 


9 





























1 


(40) 


\B\^ {E + m){E-Vof/^\C\^ |5f {E + m){E-Vof/^\C'\^ 

^ ^ (E-1/o + m)EV2 ^ + (E-1/o + m)EV2 

which agrees with the calculations of the transmission and reflection coefficient. Figure 2 shows 
the plot of the transmission and reflection coefficients. We can see from this Figure that when 
the energy of the incident electron is close to the height of the barrier the electron is completely 
reflected. As the energy of the electron increases relative to the barrier the reflection coefficient 
falls sharply and the electron is mostly transmitted as spin up. There is however a small 
probability that the transmitted electron flips its spin as well. This can be seen from the blue 
line in the upper panel of Figure 2. We show the reflection coefficients for Vq = 100 keV and 
1 MeV in Figure 3. We can see that the spin of the reflected electron depends on the height 
of the potential barrier. For Vq < rrie the reflected electron is dominantly spin up (Figure 2) 
whereas for Vq = 1 MeV > me, the electron is mostly reflected with spin down. The upper 
panel of Figure 3 shows an intermediate situation when the probabilities of reflection of electron 
as spin up and down are comparable. 


3.2 Case II: E <Vq 


We next analyze the case when the energy of the incident electron is less than the height of 
the barrier. The equations for incident and reflected electron in region I remain the same as 
equations (19) and (20). For the transmitted electron in region II, however, the equation is 
now given by 

ip 2 Ip =[—(Vo — E)iri + mir]^] ip (//) (41) 

The operator P = — (Vq — E)iri + mir]^ has eigenvalues ±iy/2(Vo — E)m = ±fp 2 - VVe choose 
the eigenvalue ip '2 which imply a decaying wave pip = = M(p)e“*'®*e“^^^) within the 

barrier. We choose the eigenvectors corresponding to eigenvalue ip 2 for spin up and down 
electrons in region II. The wave function in region II is therefore given by 
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where p = l/(Vo — E — m) and p '2 = ^/2m{E — Vq). From the continuity condition (24) we 
hnd the coefficients for this case as well. These are presented in equations (A-5), (A-6) and 
(A-7) and (A-8) of the Appendix. The reflection coefficients = |i?/Ap and R '2 = \B'/A\‘^ 
coefficients are given by 


R '2 
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(43) 

(44) 
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Figure 4: Figure shows the plot of the reflection coefficients for the case when energy of the 
incident electron is less than the height of the barrier, E < Vq. The coefficients are given in 
equations (43) and (44). For very small values of energy the incident spin up electron is reflected 
without a change in its spin. As the energy of the incident electron increases the probability 
that it flips its spin upon reflection increases considerably. 


The sum of these coefficients is always equal to one. The current densities for the transmitted 
wave (J^ and J^) is zero in this case and therefore there is essentially no wave transmitted in 
the barrier. Moreover, the reflection coefficients are independent of the height of the barrier 
for E < Vo. The probability of the reflection of spin up and down electrons is shown Figure 4. 
We can see that for very small energies the incident spin up electron is reflected with the same 
spin. As the energy of the electron increases, the probability that it is reflected with the spin 
flipped increases considerably. 

4 Extending to 3 Dimensions 

In the previous sections we saw how solutions of equation (1) leads to the transmission and 
reflection coefficients of the spin up and down electron scattering for the case of an electron 
scattering off a potential step. The analysis was performed in 1 dimension. For 3 dimensions 
the analysis is more subtle and here we discuss a version of equation (1) for 3 dimensions. 
Following is the 3 dimensional version of equation (1) 

— ifj-idiip = {irjdt + ri^m)ijj (45) 
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Here, ^idi = fiidi + ^ 2^2 + I^ 2 d 2 {i = 1, 2, 3) and the matrices /ij are given by 
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These matrices are Hermitian and obey the following anticommntation relation 


( 46 ) 

( 47 ) 

( 48 ) 


^ 2(5i,/ (49) 

Squaring the 3D equation and using (49) along with the properties of the q matrices then leads 
to the Schrodinger eqnation in 3D. The continnity eqnation can be obtained nsing a similar 
procednre employed in section 2. For this case we hrst mnltiply ' 0 ^ 73(7 + 7 ^) from the left of 
(45) and then mnltiply [t] + ?7^)73'0 to the right with the complex conjngated eqnation. Next, 
nsing the relations [hi, 73 ] = 0, {/ii ,7 + 7 ^} = 0, [ 77 ^ 73 ] = 0 and {q + 7 ^ 73 } = 0 we obtain 
the following continnity eqnation 

( 50 ) 

where Sj = + 17^)73 and T = are Hermitian matrices. The probability and cnrrent 

densities for this case are therefore given by 

Ji = (51) 

p = (52) 


Next, as before, we can seek plane wave solntions and obtain the 3D eqnation in momentnm 
space as 

fJ-iPii’ = {qE + (53) 


which is the relevant eqnation for stndying 3D scattering problems. The application of this 
eqnation to snch problems in 3D is more involved and we leave it for fntnre investigation. 

Lastly, we comment on the Panli eqnation which describes the interaction of a spin 1/2 
particle with an external electromagnetic held. It is obtained in the non-relativistic limit of the 
Dirac eqnation by assnming the presence of an electromagnetic held. Furthermore, it correctly 
predicts the spin of the particle and the gyromagnetic ratio. The Schrodinger equation, however, 
does not predict the spin of the particle. The fundamental hrst order equation introduced in 
this article is non-relativistic and it also does not predict the spin of the particle. Moreover, the 
Pauli equation cannot be obtained from this equation. However, this equation does include the 
spin of the particle in scattering problems which cannot be implemented using the Schrodinger 
equation. 
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5 Conclusion 

We proposed that the Schrodinger equation can be derived from a more fundamental hrst order 
equation in one and three dimensions, given in (1) and (45). Nilpotent symmetric matrices that 
obey anticommutation relation (5) were introduced for this purpose. We further constructed the 
solutions of this equation in ID and employed these to solve the problem of a spin up electron 
scattering from a potential step. We showed that the quantum mechanical transmission and 
reflection coefficients are obtained as the sum of those for the spin up and down electron. 
Experimental tests are required to examine the predictions of this analysis. Moreover, further 
investigation is also needed by applying this equation to other similar problems in quantum 
mechanics. 


Appendix A 

Transmission and Reflection Coefficients for E > Vq 
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Coefficients for E <Vq 
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